% 生成三角波
fs = 1000; % 采样频率
t = 0:1/fs:1; % 时间向量
f = 5; % 三角波频率，设置为低于高通滤波器的截止频率
triangle_wave = sawtooth(2 * pi * f * t, 0.5); % 生成三角波

% 设计高通滤波器
fc = 50; % 截止频率
[b, a] = butter(4, fc / (fs / 2), 'high'); % 4阶巴特沃斯高通滤波器

% 对三角波进行高通滤波
filtered_wave = filter(b, a, triangle_wave);

% 绘制结果
figure;
subplot(2, 1, 1);
plot(t, triangle_wave);
title('原始三角波');
xlabel('时间 (s)');
ylabel('幅值');
grid on;

subplot(2, 1, 2);
plot(t, filtered_wave);
title('高通滤波后的三角波');
xlabel('时间 (s)');
ylabel('幅值');
grid on;
